Weak Gravitational Lensing of CMB - Revisit 
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ABSTRACT 

In this paper, we revisit the weak gravitational lensing of CMB. The key point 
is that the deflection angle of CMB photon varies with the redshift. From this 
we obtain a set of modified power spectra. Compared with the ones calculated 
before, the temperature, the polarization E mode and the cross power spectra 
correction are of the order 0.01%, while the polarization B mode correction is 
about 0.4% (to / = 2000). 

Subject headings: cosmology: theory - cosmic microwave background - gravita- 
tional lensing 



1. Introduction 

The Cosmic Microwave Background (CMB) observation is an important probe of the 
cosmological physics, from which we have learnt much about the early universe, and obtained 
a set of precise cosmological parameters. However, in order to further constrain the cosmo- 
logical physics, we need more precise CMB measurements. In this case the secondary effects 
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g is such an effect ( 


Seliak 


1996; 


Lewis & Challinor 


2006|) (and 



In addition to the well-known Sachs- Wolfe effect and the integrated Sachs- Wolfe effect, 
the CMB photons will be distorted by the gravitational potential, after they decoupled at 
the last scattering surface. The bend of light due to weak gravitational lensing is small 
and the rms deflection is about several arcmins. But it is important at small scales, where 
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there is only little unlensed power spectrum. It is also important for polarization B mode 
( ISeljak fc Zaldarriagalll997l : iKamionkowski et al.lll997a[ ). since all the B mode power spec- 
trum comes from the weak gravitational lensing when the tensor scalar ratio is zero. This is 
because without weak gravitational lensing the scalar primord i al perturbation does not pro - 
duce any B mode power spectrum (IZaldarriaga fc Seljaklll997l : IKamionkowski et al.lll997bl ). 
However, the weak gravitational lensing will transform the polarization E mode into B mode. 
The primordial tensor perturbation also produces B mode power spectrum, whose property 
is important in constraining the inflationary scenario. Observation of such primordial tensor 
pertu rbation is regarded as a smoking gun of inflation (IMukhanov et al.lll992l : iLiddle fc Lyth 
2OOOI ). so it is very important to subtract the B mode power spectrum induced by the weak 
gravitational lensing. 

As far as we know, to account for deflection angle introduced by the weak gravitational 
lensing effect in the CMB power spectrum calculation, one use the following approximation: 
all the photons decouple instantaneously. In fact, not all photons are last scattered at 
the same time. The visibility function has two peaks which are respectively around the 
recombination era and reionization era. About (1 — e""^) ~ 10% CMB photons are last 
scattered after z ~ 10. For those photons last scattered at different times, their deflection 
angles are different, which will affect the weak gravitational lensing effect of CMB. Thus 
it affect the final CMB power spectra, especially the polarization B mode power spectrum, 
because for small tensor scalar ratio the B mode power spectrum mainly comes from the E 
mode contamination induced by the weak gravitational lensing, whereas for other types of 
power spectra the weak gravitational effect is only a secondary effect. 

This paper is organized as follows. In § |2]we provide a theoretic framework to take into 
account these factors. In § [3] we present the numerical results. Finally in § H] the discussion 



IS given. 



2. Weak gravitational lensing of CMB 

For a photon last scattered at the last scattering surface, its path of propagation is a 
curve. For example, it starts from the direction n + a{ri*, n) (77* is the conformal time of the 
recombination era), propagates along X and finally is observed in the direction n (see Fig[T]). 
Generally speaking, at different times the deflection angle is different. Besides, those photons 
which are last scattered at late time but are finally observed in the direction n also propagate 
along the same curve after they are last scattered. The ordinary method to consider the weak 
gravitational lensing effect of CMB is to set all deflection angles 0(77, n) = Q;(?7*,n). Our 
basic idea is to do our calculation by following the real line of sight X. In this paper only 
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scalar primord i al pe r turbation is discus sed and other secondary effects (jZahn et al.l 12005 
Amblard et al.l l2004l : ISantos et al.l |2003| ) such as nonlinear evolution of the gravitational 
potential are not considered. 




O 



Fig. 1. — Line of sight for CMB photons. The solid line is the real line of sight ^{f])- In the 
unlensed case it is approximately along the direction — n (Xo(?7)). 0(77*, n) is the deflection 
angle at the recombination era, while a{ri, n) is the deflection angle at a given time rj. 



2.1. Temperature 



The unlensed tempera ture anisotropies can be represented as a line of sight integral 
(ISeljak fc Zaldarriagalll996l ) 



T(n) 



_fk_ 

(27r)3 



^(k) At(?7o, 



where ^ (k) is a random variable used to describe the amplitude of the primordial fluctuation, 
which has the statistical property (^(k)^*(k')) = {27r)^P^{k)6^(k — k'). At is 
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drie-''"'STiv,k), 
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X = k{r]Q — ri), /i = k- (— n) and S'r(?7, k) is the temperature transfer function. So one obtains 

r(n) = dv j ^^a^)e-'^ST{v,k) 



VO 



dr]T{r],Xo{v)), (3) 

where Xo(?7) is the hne of hght in the unlensed case. After weak gravitational lensing effect 
is considered, the hne of sight is not Xo(?7) but instead a new curve X(r7). Thus the final 
temperature anisotropies at direction n should be 

r(n) = / dvTiv,Xiv)), (4) 



with T(?7,X(?7)) being calculated as the one in the unlensed case. It should be noted that 
generally speaking T[r), X{t))) is not the one in the unlensed case, however considering that 
the temperature source term is only low pole angular dependence, one can approximately 
replace the real T{ri,X.[rj)) with the one in the unlensed case (the same approximation is 
made in the polarization case). The new line of sight X(77) = x(n + a{ri,n)) with x = 
r]o- T] -2 drj^i, where ^ is the 'Weyl potential' defined by ^ = (^^, $^ is 

the gravitational potential appearing in the line element d'^s = a'^{ri){~{l + 2^! N)d'^ri + (1 + 
2$7v)[c?^X + fkMid"^^ + sin^ 6*6?^ </))]}). In the following the radial delays are ignored, so 
X = Vo ~ V- (^{Vj n) is the defiection angle at time i] 

a{r],n) = Vni^{r],n) (5) 

with 

V'(r/,n) = -2 /" dx ^ j^ (x'n; r/o - x'), (6) 

Jo Jk[x)Jk[x ) 



where Vn represents the angular derivative and fxix) = X ioi a. fiat universe (lKaiserlll998 



Lewis fc Challinorll2006l ). It is obvious that a{ri,n) is different from a{rj*,n). The ordinary 



method corresponds to replacing a{r], n) by a{ri*, n). 

In the flat-sky approximation, we use the 2D Fourier transformation convention 

T(r/, X(r/)) = T(r^, n + a(r^, n)) = | ^T(r^, l)e^Hfi+a(.,n))_ 
Following the standard methods, the temperature correlator between two directions can be 
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represented as 



e(r) ^ (r(n)T^(n')) 



'70 



VO 



drj I drj 



dH' 



_,il-(n+Q!(»7,n)) 



-il -{h +a(ri ,n )) 



(27r)2 J (27r)2 



VO 



VO 



drj / drj 



dH 



ilr / ^il(a(r],n)-a{r} ,n ))\^QTT{riri ) 



re (e 



where we have ignored the (T\E') term since it is very small, and in the last line we have 
already used {T{rj,\)T*{rj = {2iTf6\l-\')C'[^^'"^'^ and r = n-n. Detailed calculation of 
(jTT{r}r] ) presented i n Appendix B. (e^''("(^'")~"(^ ))) can be calculated fo llowing; the same 
methods as those of ( Lewis fc Challinor 2006 : Zaldarriaga fc Seljak IQQsh . the difference is 
that a is now time dependent. For the Gaussian variate 1- {a{rj, n) — a{rj\ n')), one has 



^i\-{a{ri,n)—a{ri ,n )) 



) = exp(- {a{rj, n) - a{rj' , n))]'^)). 



(9) 



Now one proceeds to compute the correlation tensor 



{ai{rj,n)aj{rj ,n)) = {Viilj{rj,n)Vjip{rj,n)) 

dH 



7.7 .Ji-r^i'i'ivv ) 



(10) 



where we have used {tlj{rj,\)^*{rj\\')) = {2Tr)^6\l-l')Cr^'^'^ K The detailed form of 
is given in Appendix B. One note that for a given rj and rj' , the correlation tensor can only 
depend on 6ij and the trace- free tensor f(jfj) = f^fj — j6ij. Therefore, 



{ai{rj,n)aj{7]',n')) = -CJ {r)6ij - ^^^^'^2^? (*fj) 



with 



dl 



l3cf^ivv)j^^l^^ 
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dl 



27r 



(12) 



Thus 



([1- (a(r/,n) -a(V,n))]2) 

lilj{{a{rj, n) - a{rj' , n))i{a{rj, n) - a{rj' , n))j) 



C];'{r)]+l'cos2<pC';^,{r) 



(13) 
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where we have used the following properties: 0^12(0) = 0, since ^2(0) = 0; C^^ (r) = 



V V/ 
gl \ 



Cj^gl'")) since cf^'"'^^ •* is symmetrical in the time variable and (p is the angle 



between I and r. Putting them together, one obtains the correlation function as follows, 

r-770 /-r/o , r dH 



drj I drj 
Jo 



(2vr) 



1 , c7(o) + (0) 

exp(--/2[^^i^^-^LAi _ + COS 20c;,';, (r)]) 



■no 







dr^ I dr, / :z^icj^^^^ ) 



dl 
2^' 



2"' 



exp(-^/^[ _ (,)]) 



(14) 



[•/o(/r) + -C;,»J2(/r) + ...]. 
When all deflection angles a{rj,Yv) are replaced by the value at the last scattering surface 
a{rj*, n), exp(— C^]' (r) + cos 20C^^''2(^)]) becomes a constant, which 



can be extracted out from the time integral. Thus drj dr] 
perature correlator ^(r) reduces to the ordinary one. 



TT(riri ) 
I 



CT . The tem- 



2.2. Polarization 



As was done in (IZaldarriaga fc SeljakI 119971 : iLewis fc Challinorl l2006l ) . the polarization 
field P = Q + iU {Q, U are the Stokes parameters) can be decomposed into 



+ iBk) 2Qk 

k 



(15) 



where 2Qk is the spin 2 harmonics, E^^Bk are their harmonic coefficients, and Ek,Bk are 
the harmonic coefficients of spin zero quantity E, B (corresponding to -E(n), B{h) defined in 
( IZaldarriaga fc SeljakI 119971 1 ) . These two kinds of harmonic coefficient are related by some 
quantity Sk, which expressed in spherical harmonic basis is Si = [{I — + 1)(/ + 2)]~^/^ . 
Usually E,B is the line of sight integral over time E = J^" dr] E{ri,n), B = driB{ri,n). 
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For example, scalar type perturbations produce a E fiLin fc Wandeltll2006l ) 



^ n 



vo 



drj 



d'^\< 3 



dr]E{r],Xo{v)) 



(16) 



and zero -B, where Ap is the polarization anisotropy, g{ri) is the visibility function and 
n = At2 + + Apo (At2 is the quadrupole of At, Apo(2) is the monopole (quadrupole) 
of Ap). Thus P can be rewritten as 



di] ^SkEk{ri) 2Qk 

k 

drj '^Ek{T]) 2Qk 

k 

dr]P{ri,n), 



vo 



vo 



(17) 



where Ek{ri) is the harmoni c coefficient of E(n, 'X.n(r] )) and Ek{ri) = SkEkirj). In the flat-sky 
approximation, it becomes (iLewis &: Challinorll2006l ) 



P(r7,n) 



dH 



2i(h\ „il-n 



E{r], I)e''*'e 



'1^ 



where 0i is the angle between 1 and the flxed basis e^^ in the flat-sky approximation. Analo- 
gous to P, P* can be decomposed by spin —2 harmonics. Then one obtains 



P*(r7,n) 



dH 



2i<pi il-n 



(19) 



When weak gravitational lensing effect is considered. 



P 



vo 



drj P{ri, n + a{ri, n)). 



(20) 



Now one can calculate the pol arization correlation function, as was done in the standard 
reference ( Lewis fc Challinorll2006l ). r = n— n' is chosen as the basic coordinate axis to deflne 
the Stokes parameter. The polarization fleld P^ deflned above is connected by a rotation to 
the original polarization field defined in e^^ basis, Pr = e~^*'^''P, where 0r is the rotational 
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angle. Then one obtains the correlator 

= / dr] j dr]\P{r],h + a{r],n))P{r)',n +a{r]' ,n))) 
Jo Jo 

Jo Vo V {2nyJ {2nr 



and 



g-a .(n +a{n ,n ))yi4>i^-2i^,, ^^(^^ l)E*{ri', l')) 
exp(--/ C;/ (r)]) 

[Mir) + ^^Cj^',{r) Mir) + ...], (21) 



e_(r) ^ (e-^^<^^P(n)P(n')) 

^ I dr] I dri' {e-^''^''P{ri,h+ a{ri,n)){P*{ri' ,n + a{ri' ,n)))*) 
Jo Jo 

= r dn r dn' I "^'^ / a.(n+a(^,n)) 



g-il .{n +a{v ,n ))^g-4i<^rg2i0,g2i0,, ^^^^^ \)E*{r]\ l')) 

exp(-l/2[^2^^1±^&^ - C-'(r) + cos20C-,(r)]) 

r j; dv J ^icr^-'^ 

exp(--/ ^-^ CJ/ (r)]) 



[^4(/r) + ^C^2W(^2(ir) + Mir)) + ...], (22) 



ex(r) = (e-^'^^T(n)P(n')) 

rno rvo 

= dr] rfV(e-2'*-T(r/,n + a(r7,n))(P*(r/',n +a(V,n')))*) 

Jo Jo 

vo rvo r ^2i /• j2i' 

g-a'.(n'+a(„',n')))e-2^<^-^e''*i'(-l)(T(r/,l)E*(r/',l')) 



''?o rvo 



dl 



Jo i 

exp(--/ [-^ ^ Cj;' (r)]) 

[^2(/r) + ^C;,yr)(Jo(/r) + Uh)) + ...]• (23) 

In deriving the above three formulas we have ignored the (T\l/) and {E'^) terms, and used the 
relationship = 0i-0r, (^(r], l')) = {2Tif5^{\-]^)Cf^^'^'^'^ and (T(r/, l)E*(r/', l')) = 

(27r)252(i _ 1')^^^"'' \ As in the temperature case, when all the deflection angles a{ri, n) 
are set to a(?7*, n), the three polarization correlators also reduce to the ordinary ones. 

In the above we derived the four correlators ^-(^)? from which one 

obtains the power spectra 



Cf + Cf = 21, jrdrMlr)Ur) 



Cj"" = 27r / rdrJ2{lr)U{r). (24) 



3. Numerical results 



In this paper we use the Boltzmann code CMBFAST ( ISeljak fc Zaldarriagalll996l ) and 
adopt a flat ACDM model with cosmological parameters: Qf, = 0.045, Qc = 0.225, Q/<^ = 
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0.73, Ho = 70, r = 0.08, Us = 1.0 and = 2.35 x 10"^°. When we do the integral over 

7] and 7]' , cf'^^^'^ ■* is calculated as follows: 27 points r7(j) are chosen in the region [T]st,f]o] 
('7(1) = Vst, '7(27) = Vo)y where rjst is the starting time of r] which is chosen early enough to 
ensure the contribution before it can be ignored due to the large optical depth, the same is 
done for i]'; at these 27 x 27 places (?7(j), r^^^.^), cf'^^'^^ ^ is directly calculated by Eq. (]A8l) : at 

the places where rj and rj' are around the recombination era, even more cf^^^^ ^ are directly 

calculated by Eq. (]A8|) : at other places cf^^^^ ^ is obtained by interpolation. Since cf^^^^ ■* 
vary slowly with time, the above approximation is good enough for our purpose and simplify 
the numerical calculations. 

In order to reflect the new effect, we also calculate the len sed power spectra using the 



method presented in §4.2 and §5.3.2 of (iLewis &: Challinorll2006l ) (where the flat sky approx- 



imation and correlation function method are used but all the deflection angles a{ri, n) are set 
to a;(r7-^*^, n)). Here t]'^*^ is flxed at three values r]^'^^ = 270.6Mpc, 285.4Mpc, 300.1Mpc. The 
three times are around the decoupling era and at 1]^^^ = 285.4Mpc the visibility function 
reaches its maximum. The flnal numerical results are plotted in Figs El [3], HJ [5l [61 In Figs [21 
[31 [a AC/'^-'^ is deflned as the power spectrum calculated at the deflection angle flxed at 
^fix j-Qinyg ^i^g power spectrum calculated in our paper and cj'^''^^ is the unlensed power 
spectrum. In Figs [SI [SI C^^^^ is the B mode power spectrum with the deflection angle flxed 
at Tj^^^ and Cf is the one calculated in our paper. 

The numerical results have some features. The flrst is that whatever r/-^*^' is, one cannot 
completely reproduce the real case, and generally speaking the spectra with t]^'^^ = 285.4Mpc 
are the best approximation which also deviate from the real case. For r/-^"' = 285.4Mpc, the 
B mode power spectrum change is of the order 0.1%, while other power spectra corrections 
are of the order of 0.01%. Simply speaking the effect on other types of power spectra is about 
10% of that on the B mode power spectrum. This is what is we expected because the other 
power spectra corrections due to weak gravitational lensing are of the order of 10% at around 
/ = 2000 compared with the total power spectra. The second is that the B mode power 
spectrum correction always decrease. This is because the larger is /, the more contribution 
comes from the early time. The highest relative change of B mode power spectrum minus the 
lowest one is about 0.4%. The third is that the power spectra corrections (except the B mode) 
all intersect at some / approximately. We have checked this by flxing r^-^" at some other times 
around the decoupling era and flnd that this property still holds. Besides, we also calculated 
the B mode power spectrum in the cosmological model with the same parameters but without 
reionization (in this model 77-^*^ is still flxed to be r/^^^ = 270.6Mpc, 285.4Mpc, 300.1Mpc and 
at T)-^^^ = 285.4Mpc the visibility function still reaches its maximum, see Fig 6). Comparing 
with the B mode power spectrum in the cosmological model with reionization, one can flnd 
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that the B mode power spectrum in the cosmological model without reionization do not drop 
fast at large scale, which indicates that the fast drop of B mode power spectrum at the large 
scale in the reionization case is indeed produced by reionization. 



4. Discussions 



In this paper we consider that the deflection angle of CMB photons varies with time. 
Such an effect will affect the CMB angular power spectra, especially the polarization B mode 
power spectrum, which changes about 0.4%. The smaller is the tensor scalar ratio, the more 
important is this effect. Even when the tensor scalar ratio gets as large as 0.1, the changes in 
B mode power spectrum induced by this effect cannot be ignored at small scales. If one want 
to use CMB data to provide stringent constraint on the inflationary physics, for example 
the power spectrum index of the primordial tensor fluctuation, this effect may need to be 
considered. 



We would like to thank Dr Sun Wei-Min and Qi Shi for improving the manuscript. 



A. Correlator 

In this Appendix we calculate cf^^^^ \ cf^^^^ \ cf^^^^ ^ in the case of scalar primordial 
perturbations. First we note that 



/rl^k 
t^s^(^)^t{v, k)47r J2 i'ji(x)Y:^(k)YUn). (Al) 

Im 

In deriving the last line we have made use of the formula e'"*'^ = Air Ylim ^'jK^^)^m(i^)^Zm(x). 
With the above formula, it is easy to obtain the harmonic components of T{r), n) 

/d^k 
j^^ik)STiv, k)A7rt'Mx)Y*Jk). (A2) 

After performing the angular integral, one obtains 

(Timm'm'iv)) = (A3) 

with 

Cfnvv) ^ J k^dk^P^{k)ST{v, k)ji{x)STiv',k)jiix'), (A4) 



- 12 - 



where (^(k)^*(k')) = {2'rr)^P^{k)5^(k — k') has been used. In the flat sky approximation it 
becomes 

{T{v, 1)T*(V, 0) = i2nf5\l - O^r^^^'). (A5) 

From Eq.([S]) one have 

where T^(?7i,k) is the transfer function of 'Weyl po tential'. Since \l/ = -^7 ^0 "^^ ^ Tq,(r]i,'k) 
can be directly computed by the CMBFAST code (ISeljak fc Zaldarriagalll996l ). Following 
the same steps as above, one obtains 

Wiv\ 1')) = i2nmi - (A7) 

with 



I I 



XXi Jo X Xi 

j edk^P^{k)T^{r]i,k)jiixi)T^iri[,k)ji{x[). (A8) 

From Eq. f|T6l) one have 

E{v,n) = J ^,a^)l9{v)Il{v,k){l + diy{x'e-'^). (A9) 
After doing the same thing, one can get the following formulas 

(E(r^,I)E*(V,l')) = (27r)^52(I-l')Cp('"'') 
{EME*{v,^)) = (27r)^5^(I-l>fCp('"'') 

= (27r)252(i_i')cf^(W) (AlO) 

with Si = [{I - !)/(/ + !)(/ + 2)]-^/2 and 

= Sf I kMk^-P,{k)\g{r^mr^^ fc)^-^|,(r^')n(r^', k)'^^. (All) 
For cross terms we also have 



(r(ry,l)E*(V,l')) = (27r)V(l-l>,Cp('"') 



(r(r/,l)E*(r/,l)) = (27r)V(I-l)C/ 

(27r)V(l-l>^( 

(27r)252(i_i')cf"W''^ (A12) 
with 

^Ti.(W) ^ /■ A;2rfA;-P^(A:)5T(r7, A:h(x)^^?(V)n(V, A:)^. (A13) 

J IT 4 x ^ 
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Fig. 2. — Dependence of /Cl on /. AC/ is the power spe ctrum with the deflectio n 
angle fixed at t]-^*^ (computed by the method presented in §4.2 of Jlewis fc Chalhnoilbooeh ) 
minus the one calculated in our paper, while Cf is the unlensed power spectrum. The 
solid line, long dashed line and short dashed line correspond to rj^^^ = 270.6Mpc,ri^^^ = 
285.4Mpc, T]-^*^ = SOO.lMpc respectively. 
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Fig. 3. — Dependence of AC^/Cf" on /. ACf is the power spe ctrum with the deflectio n 
angle fixed at t]'^*^ (computed by the method presented in §4.2 of Jlewis fc Chalhnoilbooeh ) 
minus the one calculated in our paper, while is the unlensed power spectrum. The 
solid line, long dashed line and short dashed line correspond to r]-^"^^ = 270.6Mpc, rj^^^ = 
285.4Mpc, r]-^*^ = SOO.lMpc respectively. 
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Fig. 4. — Dependence of l\Cj^ jiCjCf^l'^ on /. ACf^ is the power spe ctrum with the de- 
flecti on angle fixed at r]^'^^ (computed by the method presented in §4.2 of (ILewis &: Chalhnor 



20061 )) minus the one calculated in our paper, while Cj'^ is the unlensed power spectra. 
The solid line, long dashed line and short dashed line correspond to r]^'^^ = 270.6Mpc, t]^^^ = 
285.4Mpc, r]-^*"^ = SOO.lMpc respectively. 
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Fig. 5.— Dependence of (Cf ^^"^ - Cf )/Cf on /. Cf^'^ is the power spe ctrum with the de- 
flecti on angle flxed at rj-^^^ (computed by the method presented in §4.2 of (iLewis &: Chalhnor 



20061 )). while C; is the power spectrum calculated in our paper in the cosmological model 



with reionization. The solid line, long dashed line and short dashed line correspond to 
^fix ^ 270.6Mpc, r]^*^ = 285.4Mpc, r/-^^^ = SOO.lMpc respectively. 
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Fig. 6.— Dependence of (Cf ^^"^ - Cf )/Cf on /. Cf^'^ is the power spe ctrum with the de- 
flecti on angle flxed at rj-^^^ (computed by the method presented in §4.2 of (iLewis &: Chalhnor 



20061 )). while C; is the power spectrum calculated in our paper in the cosmological model 



without reionization. The solid line, long dashed line and short dashed line correspond to 
^fix ^ 270.6Mpc, r]^*^ = 285.4Mpc, r/^^^ = SOO.lMpc respectively. 



